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Abstract 

We propose and analyze the possibility of performing two-photon continuous- 
wave Doppler-cooling of hydrogen atoms using the IS — 2S transition. 

"Quenching" of the 2S level (by coupling with the 2P state) is used to increase 
the cycling frequency, and to control the equilibrium temperature. Theoreti- 
cal and numerical studies of the heating effect due to Doppler-free two-photon 
transitions evidence an increase of the temperature by a factor of two. The 
equilibrium temperature decreases with the effective (quenching dependent) 
width of the excited state and can thus be adjusted up to values close to the 
recoil temperature. 
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Laser cooling of neutral atoms has been a most active research field for many years now, 
producing a great deal of new physics. Still, the hydrogen atom, whose "simple" structure 
has lead to fundamental steps in the understanding of quantum mechanics, has not yet been 
laser-cooled. The recent experimental demonstration of the Bose-Einstein condensation of 
H adds even more interest on laser cooling of hydrogen |TJ. One of the main difficulties 
encountered in doing so is that all transitions starting from the ground state of H fall in the 
vacuum ultraviolet (VUV) range (121 nm for the IS" — 2P transition), a spectral domain 
where coherent radiation is difficult to generate. In 1993, M. Allegrini and E. Arimondo have 
suggested the laser cooling of hydrogen by two-photon tc pulses on the 15 — 3S transition 
(wavelength of 200 nm for two-photon transitions) [Q]. Since then, methods for generation 
of CW, VUV, laser radiation have considerably improved, and have been extensively used 
in metrological experiments ||. This technical progress allows one to realistically envisage 
the two-photon Doppler cooling (TPDC) of hydrogen in the continuous wave regime, in 
particular for the 15* — 25" two-photon transition. 

Laser cooling relies on the ability of the atom to perform a great number of fluorescence 
cycles in which momentum is exchanged with the radiation field. It is well known that 25* 
is a long-lived metastable state, with a lifetime approaching one second. From this point 
of view, the IS — 2S two-photon transition is not suitable for cooling. On the other hand, 
the minimum temperature achieved via Doppler cooling is proportional to the linewidth of 
the excited level involved on the process a result that will be shown to be also valid for 
TPDC. From this point of view 2S is an interesting state. 

In order to conciliate these antagonistic properties of the 1S — 2P transition, we consider 
in the present work the possibility of using the "quenching" of the 25" state to control 
the cycling frequency of the TPDC process. For the sake of simplicity, we work with a one- 
dimensional model. We write rate equations describing TPDC on the IS — 2S transition 
in presence of quenching. The quenching ratio is considered as a free parameter, allowing 
control of the equilibrium temperature. The cooling method is then in principle limited only 
by photon recoil effects. 
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We also develop analytical approaches to the problem. A Fokker-Planck equation is 
derived, describing the dynamics of the process for temperatures well above the recoil tem- 
perature T r (corresponding to the kinetic energy acquired by an atom in emitting a photon). 
A numerical analysis of the dynamics of the cooling process completes our study. 

Let us consider a hydrogen atom of mass M and velocity v parallel to the z-axis (Fig. [I]) 
interacting with two counterpropagating waves of angular frequency ujl with 2ul = ujq + 5, 
where ujq/2it = 2.5 x 10 14 Hz is the frequency corresponding to the transition IS — *> 2S, 
and also define the quantity k ~ 2k l = 2u>l/c. The shift of velocity corresponding to 
the absorption of two-photons in the same laser wave is A = hk/M = 3.1 m/s. We will 
neglect the frequency separation between 2S and 2P states (the Lamb shift - which is of 
order of 1.04 GHz) and consider that the one-photon spontaneous desexcitation from the 2P 
states also shifts the atomic velocity of A randomly in the +z or —z direction. Note that 
T r = MA 2 /fce ~ 1-2 mK for the considered transition (ks is the Boltzmann constant). We 
neglect the photo-ionization process connecting the excited states to the continuum. This 
is justified by the 1/E decreasing of the continuum density of states as a function of their 
energy E and by the fact that a monochromatic laser couples the excited levels only to a 
very small range of continuum levels. 

The atom is subjected to a controllable quenching process that couples the 2S state 
to the 2P state (linewidth T 2 p = 6.3 x 10 8 s _1 ). The adjustable quenching rate is T q . 
Four two-photon absorption process are allowed: i) absorption of two photons from the 
+2-propagating wave (named wave "+" in what follows), with a rate Ti and corresponding 
to the a total atomic velocity shift of +A; ii) absorption of two photons from the — z- 
propagating wave (wave "— "), with a rate T_i and atomic velocity shift of —A; Hi) the 
absorption of a photon in the wave "+" followed by the absorption of a photon in the wave 
"— " , with no velocity shift and iv) the absorption of a photon in the wave "— " followed by 
the absorption of a photon in the wave "+", with no velocity shift. The two latter process 
are indistinguishable, and the only relevant transition rate is that obtained by squaring the 
sum of the amplitudes of these process (called Fq). Also, these process are "Doppler-free" 
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(DF) as they are insensitive to the atomic velocity (to the first order in v/c) and do not shift 
the atomic velocity. Thus, they cannot contribute to the cooling process. As atoms excited 
by the DF process must eventually spontaneously decay to the ground state, this process 
heats the atoms. In the limit of low velocities, the transition amplitude for each of the 
four processes is the same. One thus expects the DF transitions to increase the equilibrium 
temperature by a factor of two. 

We can easily account for the presence of the quenching by introducing an effective 
linewidth of the excited level (which, due to the quenching process, is a mixing of the 2S 
and 2P levels) given by 

_ _ r„ 



2P 



<?r 2P (i) 



r 9 + r 2P 

with g = T q /(r q + r 2 p). This approximation is true as far as the quenching ratio is much 
greater than the width of the 2S state (note that this range is very large, as the width of 
the 25 state is about 10~ 8 times that of the 2P state). 
The two-photon transition rates || are given by: 

g (1 + 3MJ 2 
in_i2P 2 (5-nKVy + gy4 [ ) 

where n = { — 1,0, 1} describes, respectively, the absorption from the "— " wave, DF transi- 
tions, and the absorption from the "+" wave. I = I/I s where I s is the two-photon saturation 
intensity, 5 is the two-photon detuning divided by T 2 p, K = /cA/T 2 p = 0.26 and V = v/ A. 

The rate equations describing the evolution of the velocity distribution n(V, t) and 
n*(V,t) for, respectively, atoms in the ground and in the excited level are thus 

= _ [T _ l{V) + r + ri (V0] n(V, + y [n*(V -l) +n *(V+ 1)] (3a) 

dn*(V t) 

= r_i(V - l)n(V - 1, t) + T n(V, t) +V 1 (V+ l)n(V + l,t)- T e n*(V, t). (3b) 

at 

The deduction of the above equations is quite straightforward (cf. Fig |T]). The first term 



in the right-hand side of Eq. ( paf) describes the depopulation of the ground-state velocity 



class V by two-photon transitions, whereas the second term describes the repopulation of 
the same velocity class by spontaneous decay from the excited level. In the same way, the 
three first terms in the right-hand side of Eq. fl3b| ) describe the repopulation of the excited 
state velocity class V by two-photon transition, and the last term the depopulation of this 
velocity class by spontaneous transitions. For each term, we took into account the velocity 
shift (V — » V ± 1) associated with each transition and supposed that spontaneous emission 
is symmetric under spatial inversion. 

For moderate laser intensities, one can adiabatically eliminate the population of excited 
level. This is valid far from the saturation of the two-photon transitions and reduces the 
Eqs. (3a-3b) to one equation describing the evolution of the ground-state population: 



dn(V,t) 
dt 



, r-!(v) , Tx(v) 
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n(V,t) + 



i{r [n(V - 1, t) + n(V + 1, t)] + r_ x (V - 2)n(V -2,t) + T^V + 2)n{V + 2,t)} (4) 

Eq. (Q) is in fact a set of linear ordinary differential equations coupling the populations of 
velocity classes separated by an integer: V, V ± 1, V ± 2, • • •. This discretization exists only 
in the 1-D approach considered here, but it does not significantly affect the conclusions of 
our study, while greatly simplifying the numerical approach. 

Eqs.(|j) can be recast as dn/dt = Cn, where C is a square matrix and n is the vector 
(• • ■ n(—i, £),••• n(0, t), n(l, £),•■■). Numerically, the equilibrium distribution is obtained in 
a simple way as the eigenvector n eq of C with zero eigenvalue. In this way, the asymptotic 
temperature is obtained as : 

oo 

rp E i 2 n eq (i) 

f - = (V 2 ) = = (5) 

Tr E n eq (i) 
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Fig. shows the equilibrium distribution obtained by numerical simulation for 5 = —0.25 
and g = 1/3. The dotted curve corresponds to the distribution obtained by artificially 
suppressing DF transitions (i.e., by setting r = 0). As we pointed out earlier, the DF 
transitions lead to a heating effect. Doppler cooling is efficient mainly for atoms distributed 



on a range of g/K around the velocity V = ±\5\/K Q whereas Doppler-free transitions are 
independent of the velocity; all velocity classes are thus are affected by the heating. As a 
consequence, DF transitions induce a deformation of the velocity profile, specially for small 
values of g and 5, superimposing a sharp peak of cold atoms on a wide background of "hot" 
atoms. In what follows, all numerically calculated-temperatures are deduced form the width 
of the thin peak of cold atoms. 

Eqs. (|3a|) and (|3b|) or Eq. (§) have no exact solution. However, using some reasonable 
hypothesis, it is possible to develop analytical approaches. The most usual of these ap- 
proaches is to derive from the above equations a Fokker-Planck equation (FPE) describing 
the evolution of the velocity distribution. The derivation of the FPE for two-photon cool- 
ing follows the standard lines that can be found in the literature (see [0]). If |V| ^> 1 the 
coefficients in the resulting equation can be expanded up to second order in 1/\V\ (this is 
the so-called hypothesis of small jumps) . Moreover, if K\V\ <C \8\,g the resulting expression 
can be expanded up to the order V. The resulting FPE reads 

M dV + I + 2 ) dV 2 W 

where f = T-i(O) = To/4 and f' is the V^-derivative of T_i evaluated at V = 0. Multiplying 
this equation by V 2 and integrating over V one easily obtains: 

^ = -4f'<y 2 > + (4r + r ) (7) 

As {V 2 ) = T/T r , this equation shows that the characteristic relaxation time is (4P') -1 = 
{gT 2P I 2 5K)/{A5 2 + g 2 ). 

The equilibrium temperature is then given by 

T = 2f + r /2 = 6 2 + g 2 /4 

T r 2F K\5\ [ ' 

This results confirms that the Doppler-free two-photon transitions, corresponding to the 
contribution To/2 = 2f in Eq. (|S|) increase the equilibrium temperature (at least in the 
range of validity of the FPE) by a factor 2. This fact can also be verified from the numerical 

6 



simulations, as shown in Fig. [3|, where the dotted curve corresponds to the temperature 
obtained without DF transitions. As in one-photon Doppler-cooling, the equilibrium tem- 
perature is independent of the laser intensity (but the time need to achieve cooling obviously 
increases as the laser intensity diminishes). 

Note that the range of validity of the FPE is |V| ^> 1. It thus fails when the temperature 
approaches the recoil temperature (or, in other words, \ V\ ~ 1). Fig. f| shows the dependence 
of the equilibrium temperature as a function of the detuning for different values of parameter 
g. The minimum temperature is clearly reduced by the decreasing of g, up to values close 
to the recoil temperature T r . Moreover, the figure shows that the minimum temperature 
generally agrees with the theoretical predictions: it is governed both by the effective linewidth 
g of the excited state and by the detuning, the optimum value being 5 ~ — g/2 (in the range 
of validity of the FPE). A reasonably good agreement between numerical data and the FPE 
prediction within its range of validity is also observed. 

Let us finally note that an interesting practical possibility is to change the quenching 
parameter as the cooling process proceeds. One starts with a high value of g in order to 
rapidly cool the atoms to a few recoil velocities. Then, the quenching parameter and the 
detuning are progressively decreased, achieving temperatures of order of the recoil temper- 
ature. A detail study of the procedure optimizing the final temperature is however out of 
the scope of the present paper. 

In conclusion, we have suggested and analyzed, both analytically and numerically, the 
using of IS* — 2S two-photon transition together with the quenching of the 2S"-state to cool 
hydrogen atoms to velocities approaching the recoil limit. The quenching ratio gives an 
additional, dynamically controllable parameter. 

Laboratoire de Physique des Lasers, Atomes et Molecules (PhLAM) is UMR 8523 du 
CNRS et de l'Universite des Sciences et Technologies de Lille. Centre d'Etudes et Recherches 
Lasers et Applications (CERLA) is supported by Ministere de la Recherche, Region Nord- 
Pas de Calais and Fonds Europeen de Developpement Economique des Regions (FEDER). 
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FIG. 1. Hydrogen levels involved in the two-photon Doppler cooling in presence of quenching. 
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FIG. 2. Numerically calculated velocity distributions with 6 = —0.25 and g = 1/3. The dotted 
curve corresponds to the distribution obtained by suppressing Doppler-free transitions (cf. text). 
Typically, the distribution exhibits two structures: a broad background due to the atoms heat by 
Doppler-free transitions and a sharp peak of cold atoms. 
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FIG. 3. Dependence of the temperature (log scale) on the detuning. The full curve takes into 
account all two-photon transitions, whereas in the dotted curve the Doppler-free transitions have 
been suppressed. The plot shows that the effect of the latter is to increase the temperature by a 
factor of two, in agreement with the FPE prediction. 
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FIG. 4. Dependence of the temperature (log scale) on the detuning for three values of g: 0.9 
(full line) 0.5 (dashed line) and 0.09 (dotted line). The triangles correspond to the calculation 
based on Eq.(^|) for g = 0.5 and show the breaking of the Fokker-Planck approach at temperatures 
close to T r . The curve corresponding to g = 0.09 shows that the minimum temperature is very 
close to the recoil limit. 
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